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ATTENUATION  OF  PLANE  SHOCK  FRONTS  IN  ALUMINUM 

ABSTRACT 

The  attenuation  of  a  plane  shock  wave  by  an  overtaking  rarefaction 
has  been  studied  using  a  rotating-mirror  streak  camera  to  obtain  a  shock 
trajectory  resulting  from  the  impact  of  a  thin  striker.  Experimental 
data  are  presented  for  shocks  in  HOOF  aluminum  and  the  data  are  com¬ 
pared  with  results  obtained  from  two  sets  of  calculations.  The  shock 
attenuation  is  in  excellent  agreement  with  calculations  based  on  linear 
characteristics  along  which  (u  +  c)  is  constant.  The  attenuation  is 
measurably  less  than  that  calculated  from  Fowles '  equations,  which 
neglect  entropy  changes  at  the  shock  front. 
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TABLE  OF  SYMBOLS 
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sound  velocity  (c  =  bulk  modulus/density),  mm/ysec 

sound  velocity  of  the  undisturbed  media,  mm/ysec 

sound  velocity  immediately  behind  the  shock  front, 
mm/ysec 

striker  plate  thickness,  mm 
pressure,  kilobars 
time,  ysec 

time  required  for  the  shock  to  reach  the  rear 
surface  of  the  striker,  ysec 

time  required  for  the  shock  ttji  reach  the  rear  surface 
of  the  striker  and  the  lead  C  characteristic  to 
reach  the  striker-target  interface,  ysec 

time  required  for  the  shock  t<j>  reach  the  rear  surface 
of  the  striker  and  the  lead  C  characteristic  to 
overtake  the  shock  in  the  target,  ysec 

particle  velocity  at  the  shock  front,  mm/ysec 

particle  velocity  immediately  after  impact,  mm/ysec 

space  coordinate,  mm 

position  of  the  rear  surface  of  the  striker  when  the 
shock  front  arrives  at  the  surface,  mm 

position  of  the  striker-target  interface  when  the 
lead  C  characteristic  reaches  the  interface,  mm 

po._  tion  of  the  shock  front  in  the  target  when  the 
lead  C  characteristic  overtakes  the  shock,  mm 

( u  +  c ) ,  mm/ ysec 

(u  +  c  1  ) ,  min/ y  sec- 

constant  in  the  Murnaghan  equation,  kilobars 
(188.96  kilobars  for  aluminum) 


shock  velocity,  mm/ysec 


TABLE  OF  SYMBOLS 


j  initial  value  of  shock  velocity  in  the  target, 

mm/sec 

U,  initial  value  of  shock  velocity  in  the  striker,  mm/ysec 

o 

V  impact  velocity,  mm/ysec 

Y  constant  in  the  Murnaghan  equation,  dimensionless 

(4.266  for  aluminum) 

3 

p  density,  gm/cm 

p  density  ahead  of  the  shock  front,  gm/cm 

Mo 

a  dimensionless  parameter  equal  to  (u  +  c  -  CQ)/C0 

Note:  The  quantities  tg  and  x2  are  not  used  in  the  discussions  pre 

sented  in  this  report.  They  are  included  here  in  order  to 
make  the  remainder  of  the  symbols  agree  with  those  used  by 
Chou  and  Fowles. 
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INTRODUCTION 


The  impact  of  a  thin  striker  plate  on  a  thick  target  of  the  same 
material  produces  two  plane  shock  waves ,  one  propagating  back  into  the 
striker  and  the  other  propagating  forward  into  the  target.  The  initial 
strengths  of  the  two  shocks  are  equal  and  uniquely  determined  by  the  im¬ 
pact  velocity  and  the  Hugoniot  properties  of  the  material.  When  the 
shock  in  the  striker  reaches  the  rear  surface  of  the  plate  a  centered 
rarefaction  is  propagated  forward,  eventually  overtakes  the  forward 

moving  shock,  and  reduces  its  velocity.  Attenuation  of  plane  shock 

1* 

waves  in  solids  has  been  treated  by  G.  R.  Fowles  in  i960  and  more  recent- 

2 

ly  by  Chou,  Sidhu,  and  Zajac  .  Fowles  assumes  that  the  change  in  entropy 
across  the  shock  can  be  neglected,  an  approximation  valid  in  the  limit 
of  weak  shocks.  Chou  and  his  collaborators  have  obtained  detailed 
graphical  solutions  of  the  characteristic  equations  and  have  compared 
these  solutions  with  analytical  solutions  based  on  the  assumption  that 
the  C+  characteristics  are  linear.  For  aluminum,  the  approximate  ana¬ 
lytical  solutions  are  in  good  agreement  with  the  detailed  graphical 

solutions.  The  assumption  of  linear  characteristics  was  used  somewhat 

3 

earlier  by  Al'tshuler  et  al  to  calculate  from  experimental  data  the 
sound  velocity  behind  a  strong  shock.  Although  details  are  not  given  in 
their  paper,  they  also  performed  numerical  calculations  to  show  that  the 
error  resulting  from  the  use  of  straight  characteristic  lines  is  small. 

The  work  described  in  this  report  was  undertaken  in  an  effort  to 
provide  an  experimental  evaluation  of  the  theoretical  approximations 
described  in  the  preceding  paragraphs.  Distance-time  trajectories  for 
shocks  propagating  into  soft  aluminum  targets  were  precisely  determined 
by  optical  methods.  The  experimental  results  are  compared  with  two  sets 
of  calculations:  one  based  on  the  approximations  of  Chou  et  al  and  the 
other  based  on  the  approximations  of  Fowles. 


*  Superscript  numbers  denote  references  found  on  page  21. 
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(DISTANCE) 


FIGURE  I-  DIAGRAM  ILLUSTRATING  THE  ATTENUATION  OF  A 
PLANE  SHOCK  WAVE  BY  AN  OVERTAKING  RARE¬ 
FACTION  ORGINATING  FROM  THE  REAR  SURFACE 
OF  THE  STRIKER. 


METHODS  USED  TO  CALCULATE  SHOCK  TRAJECTORIES 

Chou's  Approximation  for  (u  +  c)  Constant. 

Figure  1  illustrates  the  shock  attenuation  problem  in  the  physical 
(x,  t)  plane.  Impact  between  the  front  surface  of  the  striker  and  the 
rear  surface  of  the  target  occurs  at  x  =  0,  t  =  0.  The  Rankine-Hugoniot 
relations  for  mass  and  momentum  conservation  together  with  the  continuity 
of  pressure  and  particle  velocity  at  the  interface  can  be  used  to  show 
that 

V 

U1  "  2  ,  (1) 

and 

Uo  =  "(Uo  “  V)  ’  (2) 

where  u  is  the  particle  velocity;  V  ,  the  impact  velocity;  U^  ,  the 
shock  velocity  in  the  striker;  and  Uq  ,  the  shock  velocity  in  the  target. 
All  velocities  are  measured  relative  to  a  stationary  observer. 

From  an  inspection  of  Figure  1,  we  note  that 
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where  and  are  the  particle  and  sound  velocities  along  the  lead 
characteristics.  When  Equations  (j)  and  (8)  are  solved  for  x^  and  t 
we  obtain 


x^  =  d 


2(Uo  ♦  cp 


2{U  -  c  ) 

0  1 


(9) 


t,  =  x_/U  (10) 

3  3  o 

after  eliminating  and  x^  by  means  of  Equations  (l)  and  (5). 

If  all  C+  characteristics  are  assumed  to  be  linear,  then  the  inter¬ 
section  of  other  characteristics  with  the  shock  trajectory  is  given  by 


x  -  x  =  z(t  -  t1)  , 


(11) 


X  - 


+ 


f  U  dt 


(12) 


where  z  =  (u  +  c)  is  a  constant  along  each  characteristic. 

For  any  given  value  of  the  parameter  z  Equations  (ll)  and  (12)  represent 
two  equations  involving  the  variables  x  and  t.  If  both  equations  are 
differentiated  with  respect  to  t  and  dx/dt  eliminated  between  the  two, 
one  obtains 


U  =  z  + 


(t  -  t1) 


dz 

dt 


(13) 


which  can  be  integrated  to  give 


t  =  t1  +  (t  -  t1)  exp 


(HO 


When  t  ,  obtained  from  Equation  (1*0,  is  substituted  into  Equation  (ll) 
one  obtains 


x 

Equations  (1*0  and 
jectory  in  which  x 


=  x  +  z  (t^  -  t1)  exp 

(15)  are  parametric  equations  for  the  shock  tra- 
and  t  are  expressed  as  functions  of  the  parameter  z. 


dz 

U-z  .  (15) 
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Fowles1  Approximation. 

In  his  analysis  Fowles  fits  the  experimental  Hugoniot  with  a  semi- 

k 

theoretical  relation  derived  by  Murnaghan  ,  which  relates  the  pressure 


and  density  by  an  equation  of  the  form 


p  =  A 


(p/p  )Y  -  1 
0 


(16) 


Fowles  used  the  Los  Alamos  shock  wave  -data^  to  evaluate  the  constants 
A  and  y.  For  aluminum  he  obtained  188.96  kilobars  for  A  and  ^.266  for  y. 
Since  entropy  changes  across  the  shock  are  neglected,  the  Hugoniot  and 
adiabat  are  identical;  hence, 


c2  =  ^i 

Po 


/  \ 
£_ 

\P°t 


Y-l 


/ 


2- 

P°l 


y-l 


(11) 


As  shown  in  Table  I,  Equation  (l6)  represents  the  pressure  data  to  with¬ 
in  1  percent  over  the  entire  pressure  range  from  a  hundred  kilobars  to 
one  megabar.  However,  the  sound  velocities  are  larger  than  the  Los 
Alamos  results  by  about  3.5  percent  at  a  hundred  kilobars  and  by  about 
1^.5  percent  at  one  megabar. 

By  combining  Equation  (l6)  with  Equation  (l)  and  making  use  of  the 
Rankine-Hugoniot  jump  conditions,  the  density  behind  the  shock  front  can 
be  related  to  the  impact  velocity. 

The  resulting  expression, 


(V/2)2  =  (A/p  )(1  -  p  /p  ) 
o  o  1 


(p ,/p  )Y  -  1 

1  o 


(18) 


is  most  easily  solved  by  graphical  or  numerical  methods.  For  t  and  x, 
Fowles  obtains  the  relations 

I - 

0 


t(o)  =  ^  +  (t3  -  t1) 


- 1 2 


Oq  -  2( y+1 ) 


-1 2 


0  -  2(y+l) 


a 


and 


x(o)  =  X  +  c  (o+l)(t  -t  ) 
1  o  3  1 


— 

0 

0 

2 

O  -2( y+1 ) 

-2(y+l) 

0 

(19) 


(20) 
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TABLE  I 


Shock  hydrodynamic  data  for  aluminum.  The  first  five 
columns  are  Los  Alamos  values  for  the  pressure,  p; 
the  relative  specific  volume,  V/V  ;  the  shock  velocity, 
U;  the  particle  velocity,  u;  and  ?he  sound  velocity, 
c,  behind  the  shock  front.  The  values  of  pressure  p 
and  sound  velocity  c  in  the  last  two  columns  were  com¬ 


puted  from 

Fowles 

equation  p  = 

188.96 

f(v  /V)1*-266  -1 
-  O  J 

p 

V 

U 

u 

c 

P 

c 

Kilobar 

V 

0 

mm/usec  mm/psec 

mm/psec  Kilobar 

mm/psec 

100 

•  9053 

6.125 

0.580 

6.307 

99-9 

6.329 

150 

.8717 

6.175 

0.831 

6.667 

150.5 

6.732 

200 

.8111 

6.793 

1.057 

6.970 

199.9 

7.092 

250 

.8211 

7.082 

1.267 

7-233 

2l9.2 

7. 123 

300 

.8007 

7-350 

1.165 

7. 165 

298.7 

7.731 

350 

.7823 

7-598 

1.651 

7.675 

319.6 

8.031 

loo 

.7662 

7.836 

1.832 

7.862 

399.5 

8.310 

1*50 

.7516 

8.062 

2.003 

8.032 

119.9 

8.576 

500 

•  7379 

8.276 

2.169 

8.190 

502.1 

8.838 

550 

•  7255 

8.180 

2.328 

8.289 

551.0 

9.086 

600 

.7113 

8.683 

2.181 

8.117 

601.9 

9.321 

650 

.70I1 

8.881 

2.628 

8.600 

655.2 

9.512 

TOO 

.6917 

9.073 

2.770 

8.713 

705.0 

9.753 

750 

.6859 

9.259 

2.908 

8.881 

751.7 

9.958 

800 

.  6779 

9.113 

3.012 

9.015 

803.2 

10.150 

850 

.6701 

9.617 

3.173 

9.  Ill 

853.5 

10.311 

900 

.0630 

9-793 

3.300 

9-271 

902.0 

10.526 

950 

.6563 

9.962 

3.121 

9-391 

950.2 

10.702 

1000 

.6198 

10.126 

3. 516 

9-508 

999-6 

10.877 

il 


where 


and 


o  =  (u  +  c  -  c  )/c 
o  o 


1  . 


(21) 


The  expressions  for  x^,  x^,  and  are  the  same  as  those  derived  in 
the  previous  section. 

Shock  trajectories  were  calculated  for  initial  conditions  corres¬ 
ponding  to  the  impact  experiments  discussed  in  the  next  section.  In  one 
set  of  calculations  the  Los  Alamos  shock  wave  data  for  aluminum  were 
used  to  evaluate  by  numerical  methods  the  integral  appearing  in 
Equations  (l4)  and  (15).  A  reproduction  of  the  Los  Alamos  data  is  pre¬ 
sented  in  the  first  five  columns  of  Table  I. 

In  the  second  set  of  calculations,  the  shock  trajectories  were  ob¬ 
tained  from  Fowles'  analytical  solution  using  the  same  impact  and  parti¬ 
cle  velocities  as  those  used  for  the  first  set  of  calculations.  The 
density  behind  the  shock  front  was  obtained  from  Equation  (l8)  and  the 
initial  shock  velocity  then  calculated  from  the  R-H  equation  for  momen¬ 
tum.  Values  of  x^,  t^,  x^,  and  t^  were  obtained  from  Equations  (5), 

(6),  (9),  and  (10)  using  values  of  c  obtained  from  Equation  (17).  The 
shock  trajectory  was  then  calculated  from  Equations  (19)  and  (20)  using 
the  value  of  obtained  from  Equation  (21). 


EXPERIMENTAL  OBSERVATIONS 

Precise  data  concerning  the  propagation  of  a  shock  can  obtained 
from  optical  records  showing  the  emergence  of  the  shock  from  the  sJant 
surface  of  a  wedge,  where  the  wedge  angle  is  chosen  sufficiently  small 
to  insure  that  the  rarefaction  from  the  slant  surface  does  not  enter 
behind  the  incident  shock  and  alter  its  velocity.  Commercially  pure 
aluminum  (HOOF)  was  chosen  for  the  present  experiments.  In  an  effort 
to  minimize  elastic-plastic  effects  associated  with  the  finite  yield 
strength  of  real  metals0,  the  aluminum  was  annealed  by  heating  it  to 
750-SOO°F,  holding  for  2  hours,  and  then  cooling  at  50°F/hour  to  500°F. 
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The  wedge  had  an  angle  of  15  degrees  and  a  square  base  whose  sides 
were  5-1 A  inches.  The  striker  plate  was  mounted  on  a  7-1/2  inch  di¬ 
ameter  charge  and  positioned  3A  inch  from  the  base  of  the  wedge.  The 
base  charge  was  initiated  by  a  6  inch  diameter  planewave  lens.  Because 
of  difficulty  in  maintaining  the  integrity  of  explosively  accelerated 
plates,  the  strikers  were  made  from  2021-T3  aluminum  alloy.  However, 
errors  due  to  an  elastic-plastic  effect  in  the  striker  should  be  negli¬ 
gible  when  the  shock  trajectory  is  measured  over  a  distance  correspond¬ 
ing  to  more  than  twenty  thicknesses  of  the  striker.  A  diagram  illus¬ 
trating  the  essential  features  of  the  experimental  arrangement  is  pre¬ 
sented  in  Figure  2. 

The  explosive  charge  assemblies  were  varied  to  provide  the  three 
combinations  of  striker  velocities  and  thicknesses  shown  in  the  table 
below: 


Striker  Thickness 


Velocity 


1/16  in.  (1.588  mm) 
3/32  in.  (2.381  mm) 
1/8  in.  (3.175  mm) 


5.800  mm/ysec 
5.130  mm/ysec 
1.136  mm/ysec 


Two  successful  firings  were  obtained  with  3/32  inch  and  1/8  inch  strik¬ 
ers.  Only  one  successful  firing  was  obtained  with  l/l6  inch  strikers 
due  to  a  malfunction  of  the  initiation  system. 


The  most  accurate  determination  of  impact  velocity  was  obtained  by 
using  the  linear  part  of  the  shock  trajectory  (that  portion  between 
x  =  0  and  x  =  x^)  to  determine  the  initial  shock  velocities.  The 
Hugoniot  data  for  aluminum  were  then  used  to  establish  the  corresponding 
values  of  u^,  which,  in  accordance  with  Equation  (l),  are  one-half  the 
impact  velocities.  The  complete  shock  trajectory  for  the  impact  by  the 
3/32  inch  plate  is  presented  in  Figure  3.  Reproducibility  of  the  shock 
trajectory  is  illustrated  in  Figure  3  by  the  excellent  agreement  between 
data  obtained  from  two  separate  impacts. 
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FIGURE  2.  DIAGRAM  ILLUSTRATING  THE  USE  OF  A  ROT  ATI  NG  —  MIRROR 
CAMERA  TO  MEASURE  SHOCK  ATTENUATION  IN  A  METAL  WEDGE. 


(  D3S rh  3WI1 


FIGURE  3- GRAPH  SHOWING  THE  SHOCK  TRAJECTORY  RESULTING  FROM  A  3/32  IN. 

STRIKER  IMPACTING  AT  5.  1 3  O  MM/ ft  S  EC.  DATA  FROM  TWO  SEPARATE  IMPACTS 
ARE  DISTINGUISHED  BY  USING  CIRCLES  FOR  ONE  AND  TRIANGLES  FOR  THE 
OTHER. 


DISCUSSION  OF  RESULTS 


Arrival  times  for  the  shock  at  various  depths  in  the  target  are 
presented  in  Table  II.  Column  A  represents  experimental  data;  values 
in  Column  B  were  obtained  from  the  Los  Alamos  sound  velocity  data  using 
the  (u  +  c)  =  constant  approximation  of  Chou  et.  al.;  and  values  in 
Column  C  were  calculated  from  Fowles  equations.  Times  calculated  using 
the  Los  Alamos  sound  velocity  data  are  in  excellent  agreement  with  the 
experimental  results,  the  maximum  discrepancy  being  only  0.08  usee.  Al¬ 
though  the  differences  are  exceedingly  small,  results  obtained  from 
Fowles'  equation  show  that  neglect  of  entropy  changes  produces  a  shock 
attenuation  which  is  more  rapid  than  that  observed  experimentally.  As 
would  be  expected  entropy  changes  become  less  important  as  the  initial 
shock  strength  decreases. 
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TABLE  II 


Comparison  of  measured  shock  trajectories  with  those 
comj  ted  using  two  different  approximations.  Column  A 
contu.ns  values  interpolated  from  the  experimental  data, 
Column  B  contains  values  of  t  computed  from  the  Los  Alamos 
sound  velocity  data  using  the  method  of  Chou,  et  al., 
and  Column  C  contains  values  of  t  computed  from  Fowles ' 
approximation  neglecting  the  change  in  entropy  across 
the  shock. 


1/16  in. 


X 

A 

B 

C 

A 

mm 

psec 

psec 

psec 

psec 

0 

0 

0 

5 

0.55 

0.57 

10 

1.10 

1.08 

1.08 

1.13 

15 

1.67 

1.70 

20 

2.29 
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